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The framework of non-relativistic effective field theory (NREFT) aims to generalise the standard 
analysis of direct detection experiments in terms of spin-dependent (SD) and spin-independent (SI) 
interactions. We show that a number of NREFT operators lead to distinctive new directional 
signatures, such as prominent ring-like features in the directional recoil rate, even for relatively low 
mass WIMPs. We discuss these signatures and how they could affect the interpretation of future 
results from directional detectors. We demonstrate that considering a range of possible operators 
introduces a factor of 2 uncertainty in the number of events required to conhrm the median recoil 
direction of the signal. Furthermore, using directional detection, it is possible to distinguish the more 
general NREFT interactions from the standard SI/SD interactions at the 2cr level with 0(100 — 500) 
events. In particular, we demonstrate that for certain NREFT operators, directional sensitivity 
provides the only method of distinguishing them from these standard operators, highlighting the 
importance of directional detectors in probing the particle physics of dark matter. 


I. INTRODUCTION 

The detection of dark matter (DM) in a laboratory 
setting is considered one of the greatest goals of mod¬ 
ern particle physics. Direct detection experiments mm 
aim to measure the keV-scale recoil energy imparted on 
detector nuclei by interactions with Weakly Interacting 
Massive Particles (WIMPs) in the Galactic halo. The 
motion of the Earth and Sun in the Galactic rest frame 
induces a highly directional flow of DM particles in the 
lab frame, known as the WIMP wind’. The result is that 
the mean recoil direction is opposite that of the Solar mo¬ 
tion, in the direction of the constellation Vela. Detection 
of this directional signal would provide strong evidence 
for a DM origin of the signal |3H5]. 

A number of experiments with directional sensitivity 
are currently in development. The most advanced of 
these are low pressure time projection chambers (TPCs) 
such as DRIFT PIT], MIMAG PP, DMTPG [13E], 
NEWAGE [121 [E] and D3 [T4| [15] , though a number of 
other technologies are also being considered, including 
nuclear emulsions m, DNA-based techniques m and 
the possibility of exploiting columnar recombination in 
Xenon targets [T11[T9]. The analysis of data from such ex¬ 
periments (as well as from non-directional experiments) 
typically assumes standard spin-dependent (SD) or spin- 
independent (SI) contact interactions [501 E] which are 
leading (zeroth) order in the recoil momentum q and DM- 
nucleus relative velocity v. This is because WIMPs in 
the Galactic halo are highly non-relativistic {v ~ 10“^c), 
leading to typical momentum transfers in the 100 MeV 
range. 

However, in recent years, a great deal of effort has 
been focused on developing a more general framework 
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for DM-nucleus interactions. The non-relativistic effec¬ 
tive field theory (NREFT; introduced by Fan et al. [55] 
and extended in Refs. [53H55] ) considers all possible non- 
relativistic quantum mechanical operators which can con¬ 
tribute to elastic DM-nucleus scattering, including those 
which are higher order in q and v. Such a basis of oper¬ 
ators allows us to avoid biased reconstructions of future 
DM signals by accounting for all possible DM-nucleon 
interactions. In addition, it enables us explore possible 
blind spots in the sensitivity of current experiments. To 
this end, a number of authors have presented constraints 
on these operators from current and future direct detec¬ 
tion [5SH5T] and neutrino telescope data [55Ij55] - 

In this work, we consider the impact of NREFT op¬ 
erators on the recoil spectra in directional direct detec¬ 
tion experiments. As we discuss in detail in Sec. jTTTJ the 
event rates arising from NREFT operators typically ap¬ 
pear with additional powers of q^ and/or when com¬ 
pared with the standard SI/SD operators, where v± is 
the DM velocity perpendicular to q. Additional powers 
of q^ typically result in a spectrum which grows with re¬ 
coil energy En. When considering the directional rate 
integrated above a certain threshold energy Emin, this 
enhances the contribution of directional features at high 
energy, leading to a more sharply directional event rate. 
In contrast, powers of tend to suppress recoils in the 
forward direction, leading to a more isotropic distribu¬ 
tion. 

These differences can have significant consequences for 
interpreting future data from directional data. Calcu¬ 
lations of the number of WIMP signal events required 
to discriminate from an isotropic background or to con¬ 
firm the median recoil direction are performed assuming 
standard SI/SD interactions. As we will show, more gen¬ 
eral directional rates may require more or fewer events to 
confirm the WIMP origin of a signal. Additionally, some 
NREFT operators lead to distinctive directional signa¬ 
tures which may allow them to be distinguished from the 
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standard scenario using directional detectors. 

The key result of this paper is presented in Fig. 
which shows the statistical significance with which stan¬ 
dard interactions can be excluded as a function of the 
number of signal events. We show that certain NREFT 
operators (in particular O 7 and Og,, see Eq. A2) can be 
distinguished from the standard SI/SD operators at the 
95% confidence level with several hundred signal events 
in a directionally-sensitive Fluorine-based detector. For 
these operators, discrimination would be almost impos¬ 
sible using only information about event energies, high¬ 
lighting the important role directional detectors may play 
in exploring the particle physics of DM. 

In Sec. |nj we describe in more detail the directional 
event rate. We then introduce the NREFT operators 
in Sec. |III| and demonstrate how their directional spec¬ 
trum is expected to differ from the standard case. In 
Section |IV) we consider a variety of directional statistics 
which have been proposed to distinguish a directional 
DM signal from isotropic backgrounds. Using these, we 
determine how many signal events are required to de¬ 
tect anisotropy and confirm the median recoil direction 
of the signal, and compare these results with the SI 
and SD rates which are typically considered. Finally, 
we discuss direct comparisons between the standard di¬ 
rectional rate and these more general interactions. We 
show how directional detection can be used to distinguish 
non-relativistic operators which couple to the transverse 
velocity-squared from the canonical SI and SD operators. 


II. DIRECTIONAL RATE 


We begin with a summary of the standard event rate in 
directional detectors. The double-differential recoil spec¬ 
trum per unit detector mass for DM-nucleus scattering 
with a fixed DM velocity v is given by [35] : 
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Here, and mjv are the DM and nuclear masses re¬ 
spectively, po is the local DM density and q is the direc¬ 
tion of the recoiling nucleus. The matrix element-squared 
(|AI| ), summed and averaged over the initial and final 
spins, is determined by the particle physics operators 
which mediate the interaction. Finally, the d-function 
imposes the kinematic constraint on the elastic scatter¬ 
ing, where r’min is the minimum WIMP speed required 
to excite a nuclear recoil of energy Epi, 
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with the DM-nucleus reduced mass given by = 

/{m^ + rriRi). 

In the standard analysis framework, the matrix ele¬ 
ments are calculated assuming contact interactions which 


are leading order in the momentum exchange and rela¬ 
tive DM-nucleus velocity. This is because virialised DM 
particles in the Galactic halo are expected to have highly 
non-relativistic speeds, v ^ 10 “^c, so any higher order 
interactions will be suppressed by factors of 10“®. These 
leading order interactions are assumed to arise from a 
coupling of the spins of the DM and nucleons (spin- 
dependent, SD) or a coupling of their particle densities 
(spin-independent, SI): 


OsD = ■ Sn 
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where we use the subscript n for nucleon. 

The resulting matrix elements are then written as 

(\M\^) = (cPQP + c^O"), (4) 

where and c" are the coupling strengths of the oper¬ 
ators with protons and neutrons. The term in angular 
brackets on the right hand side is referred to as the nu¬ 
clear response function, and is the expectation value of 
the operators (either SI or SD) over all nucleons in the 
nucleus. For the SI case, this gives a coherent enhance¬ 
ment of the scattering rate roughly proportional to 
for a nucleus composed of A nucleons. In the SD case, 
the response function takes into account the total nu¬ 
clear spin, as well as the expectation values of the proton 
and neutron spin within the nucleus. For both types of 
interaction, the finite size of the nucleus leads to a loss 
of coherence at large momentum transfer, meaning that 
the nuclear response functions give a roughly exponential 
suppression of the rate with recoil energy |37| |38| . 

The final component for calculating the standard di¬ 
rectional event rate arises from the fact that DM in the 
Galactic halo has a distribution of velocities f{v). Thus, 
we must integrate Eq. over all DM velocities, weighted 
by f{v): 
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All dependence on the velocity distribution has been ab¬ 
sorbed into the Radon Transform (RT) [33], defined as: 

/(Wmin,g)=/ f{v)S {v ■ q- Vmin) d^V, (6) 
7 r 3 


where we have changed variables in the argument of the 
(5-function, leading to an extra power of v in the inte¬ 
gral. Physically, the RT is obtained by integrating over 
all velocities for which the observed recoil is kinemati¬ 
cally allowed. 

In the Standard Halo Model (SHM), dark matter is 
assumed to follow a Maxwell-Boltzmann velocity distri¬ 
bution, given by 
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For an isotropic, isothermal sphere of DM, with den¬ 
sity profile p oc the average velocity of the DM 
particles with respect to the Earth Flag is related to 
the velocity dispersion by uiag = y/^cTy. A value of 
uiag « 220 km s“^ is typically used [JD], though values 
in the range 180 — 270 km s“^ have been suggested [4H - 
l45] . The corresponding RT is given by [36] : 
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However, we note briefly that the SHM is unlikely to 
be a realistic description of the true DM distribution. Re¬ 
sults from N-body simulations indicate deviations from 
the standard Maxwell-Boltzmann distribution |46l Wf\ . 
including the possibility of additional structures such as 
dark disks |i5H5T] or streams |S2J[^. For concreteness, 
we assume the SHM as a standard benchmark, with fixed 
values of viag = 220 km s~^ and = 156 km s~^ in this 
study. We will also neglect effects due to the finite Galac¬ 
tic escape speed |S 1 ], which will not be significant over 
the range of recoil energies considered here. We leave an 
exploration of the impact of astrophysical uncertainties 
to future work. 

It will sometimes be necessary to distinguish the full 
double-differential recoil rate d^R/dAi^dllg from the 
energy-integrated recoil rate, given by 
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In this case, we are interested in the direction of all 
recoils observed in the detectors (in an energy window 
E e [E’min, Ajiiax]), but not the energy of each event. 


III. NON-RELATIVISTIC EFFECTIVE FIELD 
THEORY (NREFT) 


In non-relativistic effective field theory (NREFT), the 
standard set of SI and SD operators are extended to in¬ 
clude all those interactions constructed from Galilean, 
Hermitian and time-reversal invariant operators |23| . 
This framework was extended in Ref. |21| to include com¬ 
posite operators which do not typically arise due to the 
exchange of spin-0 or spin-1 mediators. The NREFT in¬ 
teraction operators are rotational-invariants constructed 
from the following Hermitian operators: 
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Here, is the momentum transferred to the interacting 
nucleon, are the WIMP and nucleon spin operators 
and the operator v:^ is defined as 
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This is the component of the DM velocity v perpendicular 
to the recoil momentum. By energy conservation, we 
therefore have • ()„ = 0. The DM velocity F does 
not have definite parity under the exchange of incoming 
and outgoing particles, and is therefore not Hermitian. 
The transverse velocity v:^, however, is Hermitian and is 
therefore the only combination in which the DM velocity 
may appear. 

The full list of possible WIMP-nucleon operators is 
given in Appendix (Eq. |A2| ) , with notation matching 
that given in Refs. [231124j. Within this framework, we 
write the full interaction Lagrangian as 

C = ^c°0° + clO}, ( 12 ) 

i=l 

where the superscript indices 0,1 denote the isoscalar and 
isovector couplings and operators respectively!^ Within 
the NREFT framework, the standard SD and SI interac¬ 
tions are labeled O 4 and Oi respectively, and have the 
exact form given in Eq. [^ 

The matrix element-squared is then written as 

(1-^1')= E E clcj'El/(vl,q^), (13) 

i,j=l 

where E^ are the nuclear response functions associated 
with the and j**' operator. For simplicity, we will ne¬ 
glect interference terms (i 7 ^ j). In addition, we neglect 
the operator O 2 , which does not typically appear at lead¬ 
ing order in the non-relativistic reduction of a relativistic 
interaction Lagrangian. The list of nuclear response func¬ 
tions in terms of a set of standard form factors is given in 
Eq. |A4| of Appendix |A| We note that the standard SI and 
SD form factors are often normalised to unity at g = 0 
m- However, in this framework, we include any coher¬ 
ent enhancement factors in the definition of the nuclear 
form factors. 

We will also consider an example of a long-range DM- 
nucelon operator, as described in Refs. [S5H57j. Such op¬ 
erators are not contact operators but instead arise from 
the exchange of light mediator particles and are therefore 
not strictly effective field theory operators. However, we 
include an example in this work due to the novel de¬ 
pendence which they give rise to. The example we con¬ 
sider is = 0 \lq^ which behaves as 0 \ (the standard 
SI operator) with an additional 1/g'* suppression of the 
cross section. 

We note that the response functions Eij depend on the 
incoming and outgoing momenta only through and 
[56] ■ which can be written as 
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^ In the proton-neutron basis, the couplings can be written c® = 
-h c^) and = |(c^ — c^). 


( 11 ) 
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These are now the 'WIMP-nucleus momentum transfer 
and transverse velocity respectively. The operators form¬ 
ing the basis of the NREFT are rotationally invariant, so 
there is no preferred basis in which to measure the di¬ 
rections of v±_ and q. They can therefore only appear as 
scalar products. However, by construction ?7j_ • g = 0. 
Moreover, averaging over nuclear spins means that 
does not pick out a particular direction. The result is 
that the response functions in Eq. [T^ depend only on the 
magnitude of the recoil momentum g, but not its direc¬ 
tion g. Substituting these response functions into Eq. 
has no impact on the directionality of the signal; the only 
dependence on the recoil direction is through v ■ q within 
the d-function. 

However, in passing from Eq. to Eq. care must 
be taken, because now the matrix elements can depend 
on the DM velocity, via their dependence on v^. The 
full directional event rate, integrating over the WIMP 
velocity distribution, is then 


A. Transverse Radon Transform 


In order to calculate the TRT, we can decompose the 
DM velocity into components perpendicular and paral¬ 
lel to q, V = (17l,'C||), where we note that v± is a two- 
dimensional vector. In this basis, we can write v-q = v\\, 
meaning that Eq. [T^ reduces to 


f'^{Vnun,q) = C {v±)'^f{v±,V\\=Vmin)d^V±- ( 17 ) 


The requirement that wn = timin ensures that the kine¬ 
matic constraints of the elastic scattering are satisfied. It 
is then necessary to integrate over all possible transverse 
velocities, weighted by the transverse velocity squared. 
Geometrically, we must integrate f(v) over a plane per¬ 
pendicular to g, at a distance Wmin from the origin, 
weighted by the square of the perpendicular distance 
along the plane. For the SHM, given by Eq. the TRT 
becomes 
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The response functions Fij are composed of terms which 
are proportional either to or In the former case, 
the integral over the velocity distribution is simply the 
Radon Transform of Eq. In the latter case, we must 
compute the transverse Radon Transform (TRT), which 
we define as 


f'^{Vmin,q) = C {v±ff{v)6{v-q-Vmin)d^V. 


(16) 


We have defined the TRT as carrying two inverse factors 
of the speed of light. In standard natural units, speeds 
are dimensionless, so this ordinarily does not need to be 
made explicit. However, we include this factor in order 
for the TRT to have the same units (of inverse speed) as 
the standard RT, to allow a more transparent compari¬ 
son. 

We note that for some operators (e.g. Oj), the relevant 
nuclear response function is proportional to , meaning 
that the directional rate is proportional to the TRT of 
the velocity distribution. In other cases (e.g. Oio), the 
form factor has no dependence on and therefore the 
directional recoil rate (Eq. [T5| has the same directional 
dependence as the standard SI/SD scenario. In general, 
however, both types of terms may be present and the 
full directional rate may be somewhere between the two 
regimes. 


where = v\g,g — viag ■ q. Performing the integral over 
transverse velocities, we obtain 
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For more extreme distributions, such as a stream 
[53| . the velocity distribution can typically be modelled 
as a Maxwell Boltzmann distribution |58j with a small 
velocity dispersion and ?7iag = Vs, the stream velocity. 
However, in the extreme case that ay —>■ 0, the stream 
distribution becomes 


f{v) = {v-Vs) . (20) 


The corresponding TRT in this case is then given by 

- Vs ' q) . (21) 

In order to calculate the total energy spectrum of 
events, it is necessary to integrate the RT or TRT over 
all angles, depending on the relevant operator, in order 
to obtain the corresponding velocity integrals: 
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FIG. 1. Comparison of standard and transverse Radon Transforms for the SHM. The left panel shows the Radon 
Transform (RT) of the SHM (defined in Eq. [^, integrated over the azimuthal angle cj> and with -uiag aligned along 8 — 0. The 
right panel shows the corresponding transverse Radon Transform (TRT) which appears in the directional rate for NREFT 
operators coupling to v± and is defined in Eq. |16| 
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FIG. 2. Comparison of standard and transverse Radon Transforms for a stream. The left panel shows the Radon 
Transform (RT) of a stream distribution function, integrated over the azimuthal angle (f) a-nd with uiag aligned along 9 — 0. 
We approximate the stream as a Maxwell-Boltzmann distribution with uiag = 400 km s“^ and — 20 km s“^. The right 
panel shows the corresponding transverse Radon Transform (TRT) which appears in the directional rate for NREFT operators 
coupling to v± and is defined in Eq. |16| 


We have verified explicitly that the velocity integral 
obtained in this way matches that given in Ap¬ 
pendix A of Ref. [SS] . 

In Fig.[T] we compare the standard and transverse RTs 
for the SHM. For ease of presentation, we have integrated 
over the azimuthal direction (j) and chosen the angular 
basis such that Flag is aligned along 9 = 0. The angle 9 
is therefore the angle between Flag and the recoil direc¬ 
tion. The value of the TRT is approximately 10® times 
smaller than the RT, as the TRT is suppressed by two 
powers of u/c ~ 10“®. The most striking feature of the 
TRT, however, is that the maximum in the differential 
recoil rate does not occur along 0 = 0 as in the standard 
case. This is because forward-going WIMPs cannot in¬ 


duce forward-going nuclear recoils, due to the weighting 
by v]_. Instead, the maximum rate occurs approximately 
perpendicular to Flag, where is maximised. However, 
the finite width of the SHM distribution means that a sig¬ 
nificant population of WIMPs will have velocities which 
deviate from Flag- This means that recoils along the di¬ 
rection of Flag are still possible, so the forward scattering 
rate is not precisely zero. 

For comparison, we show in Fig. the more ex¬ 
treme example of a stream, modelled as a Maxwell- 
Boltzmann distribution with uiag = 400 km s~^ and 
ay = 20 km s~^. As in the case of the SHM, the peak re¬ 
coil direction deviates substantially from the forward di¬ 
rection. However, in the case of the stream, recoils in the 
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forward direction are almost entirely suppressed. This 
is because of the very narrow velocity dispersion, which 
means that all particles are travelling with velocity close 
to zTiag and so cannot induce recoils in that direction. 
The result is that, compared with the standard RT, the 
TRT is truncated at large values of Umin, above around 
350 km s“^. 


B. Comparing NREFT operators 



We now compare the directional rate obtained for the 
different NREFT operators. For concreteness, we will 
consider a CF 4 target, which is used in several gaseous 
TPC experiments [71 1 mi ng, and which provides a 
promising SD WIMP-proton target. We will focus only 
on interactions with Fluorine. Carbon makes up only 
14% of CF 4 by mass and is spin-zero and we therefore 
expect the contribution from Carbon to be subdomi¬ 
nant. We consider a WIMP of mass rriy^ = 100 GeV 
with SHM velocity distribution and an experimental sen¬ 
sitivity in the energy window G [20, 50] keV. An en¬ 
ergy threshold of 20 keV has previously been reported 
by the DRIFT-IId experiment m, although the angular 
resolution of directional experiments worsens at low en¬ 
ergies jSD]. We limit our analysis to isoscalar couplings 
(c^ = c"), though as we will see, the differences in direc¬ 
tionality arise predominant from the scaling of the dif¬ 
ferent response functions with u_l and q, so we do not 
expect the results to change substantially for more gen¬ 
eral couplings. 

Figure [^ shows the total directional rate integrated 
over the energy window of the experiment (defined in 
Eq.§ expressed as a function of 9, the angle between 
Flag and the nuclear recoil direction. The directional rate 
has been normalised to unity to allow a comparison of the 
angular distribution of events between different NREFT 
operators. 

We show results only for a selection of operators. We 
find that each of the remaining operators leads to a direc¬ 
tional rate which is almost indistinguishable from one of 
those plotted in Fig. For light nuclei such as Fluorine, 
form factors do not decay as rapidly with q as for heavier 
nuclei (and can often be assumed to be approximately 
constant m), so we expect that differences in form fac¬ 
tors for the different operators should not be significant. 
This grouping of different operators therefore arises due 
to their different functional dependence on powers of q^ 
and We classify the nuclear response functions for 
Fluorine as follows (in a similar fashion to the classifica- 


FIG. 3. Directional event rates, normalised to unity, 
for several NREFT operators: O3 (red A), O4 (black •), 
O5 (cyan *), O7 (blue T), Og (magenta ♦), O15 (green ■) 
and Oi^ (yellow ►). The form of these operators is given 


in Eq. |A2[ Each of the remaining NREFT operators gives 
a directional distribution similar to one of those shown here, 
depending on the functional dependence of the operator in 
question (see Eq. 241. The angle 9 is defined as the angle be¬ 
tween Flag and the direction of the nuclear recoil. We assume 
= 100 GeV and a Fluorine detector with rate integrated 
over Er G [20, 50] keV. 


tion of Ref. [24]]l: 


Proportional to 


I 

Oi,04, 

vl 

07,Os, 

g2 

Og, Oio 

vW 

05,0l3 


03,^6, 

g"‘(g^ + u5_) 

^15 , 

9 -^ 

C>f«. 


C’li, O12 , 
Ol4, 


(24) 

Operators belonging to the same class will lead to ap¬ 
proximately the same directional rate, so we therefore 
show only a single example from each class in Fig. [^ 
The standard directional signal arising from SD (or 
SI) interactions corresponds to the operator O 4 in Fig. 
This standard signal lies intermediate between the re¬ 
maining NREFT operators. Nuclear response functions 
suppressed by positive powers of lead to a directional 
rate more sharply peaked towards Flag. This may be 
surprising, given that such response functions will affect 
only the energy dependence of Eq. |I5| and not the angular 
dependence (for a fixed value of Er). However, the dis¬ 
tribution of recoils becomes increasingly anisotropic with 
increasing Umin, which can be seen from the left panel of 
Fig. El When we integrate over all energies, response 
functions which scale as give a greater weight to this 
more anisotropic recoil distribution at high Umin, leading 
to a more directionally peaked distribution. Equivalently, 
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we note that by energy conversation, q = 2 /r^jvu • q = 
cos 9. This means that increasing the relative con¬ 
tribution of recoils at large q is equivalent to increasing 
the contribution of recoils with small values of 9, leading 
to a more peaked directional spectrum. Further pow¬ 
ers of 9 ^ leads to increasingly forward-peaked directional 
spectra. 

The long-range operator has a nuclear response 
function which scales as Again, the directional de¬ 
pendence of the full double-differential rate is unchanged 
compared to the standard case. However, integrating 
over all energies, the greatest contribution now comes 
from low energy recoils, for which the spectrum is most 
isotropic. This leads to a directional spectrum which is 
less peaked in the forward direction (relative to the stan¬ 
dard SI/SD interactions), as shown in Fig. 

Nuclear response functions which are suppressed by 
powers of also lead to directional spectra which are 
less sharply peaked than the standard case. However, 
in this case, it is due to a fundamental difference in the 
directionality of the double-differential recoil spectrum. 
This behaviour is encoded in the TRT and arises because 
scattering in the forward direction is suppressed by the 
coupling to tij_, while scattering perpendicular to Flag is 
enhanced. In particular, for the operators 0^ and O 5 
(and operators with similar functional forms) the peak in 
the recoil distribution occurs at some non-zero angle from 
the direction of peak flux Fag. We can see this clearly 
in the right panel of Fig. where the TRT increases 
moving away from the forward recoil direction. For the 
operator O 5 , an additional suppression by q^ means that 
the directional spectrum peaks at a slightly lower angle 
than in the case of O 7 . However, the result for both 
operators is that the recoil spectrum will exhibit a ring¬ 
like feature, with the peak recoil direction being greatest 
in a ring around the median recoil direction. 

We illustrate the features of this ring in Fig. |^for the 
operator O 7 , as a function of the WIMP mass and thresh¬ 
old energy. We continue to consider the directional spec¬ 
trum integrated over recoil energies, defined in Eq. 
The solid contours indicate the ring opening angle in de¬ 
grees (that is, the angle between Fag and the peak recoil 
direction). The shaded regions show the ring amplitude: 
the ratio between dR/d9 at the peak and at 9 = 0. 

A similar ring-like feature in the standard directional 
rate has been studied previously ED- When considering 
the standard RT, a ring is only possible if Vmin < t'lag, 
meaning that the term inside the exponential in Eq. 
can be set to zero for some value of 9. For the TRT, 
there is also a dependence on 9 outside the exponential. 
By differentiating Eq. [T^ with respect to 9, it can be 
shown that a maximum in for 9 0 therefore exists 

for larger values of z;i„in than in the standard case, up 
to Wmin = 2wiag. This means that the ring arising from 
NREFT operators which couple to v]_ can be observed 
for lower WIMP masses, higher energy thresholds and for 
smaller values of flag- For example, for the range of pa¬ 
rameter values displayed in Fig.|^ no significant ring-like 
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FIG. 4. Properties of the ring in the directional re¬ 
coil spectrnm (integrated over Eh € [Amin, 50] keV) for the 
operator Oy, assuming a Fluorine target. The ring opening 
angle in degrees is indicated by solid contours, while the ring 
amplitnde (ratio of the maximnm rate to the rate at 0 = 0) is 
shown by the blue shaded regions. The maximum ring open¬ 
ing angle for the parameters considered here is approximately 
52°. 


feature is observable in the directional spectrum (inte¬ 
grated over energy) for the standard SI or SD operators, 
giving a maximum ring opening angle of ^ 1 ° (compared 
to 52° for the O 7 operator)]^ A detailed study of ring-like 
features is conducted in Ref. ED, but we note here that 
if a prominent ring is observed for low mass WIMPs or 
for a relatively high threshold, this may be indicative of 
non-standard operators, which couple to the transverse 
WIMP speed. 

In general then, the basis of NREFT operators leads 
to a variety of directional signatures. Those operators 
suppressed by powers of q'^ lead to a more anisotropic 
directional spectrum (when we consider event directions, 
integrating over all event energies), while those which 
couple to give a more isotropic spectrum, with the 
possibility of a ring-like feature in the spectrum, even for 
relatively large values of Umin- 

In the rest of this work, we focus on the operators O 7 
and O 15 , which are written. 


Or = S„-v^, (25) 

Oi 5 = -(4 • —)((4 X Fi) • ^), (26) 

and lead to the following nuclear response functions, 

Ft ,7 cx vlFs ^, (27) 


^ A significant ring is found in Ref. m for larger values of 
than those assumed here or, alternatively, when the full spectrum 
is considered (i.e. when we do not integrate over energy). 










F^„j . (28) 

As can be seen in Fig. the operator O 7 leads to the 
most isotropic directional rate, while the operator O 15 
leads to the most sharply peaked directional rate. We 
note in the case of O 15 that there are two terms in the 
nuclear response function, one proportional to q® and the 
other proportional to q'^v^. These give similar contribu¬ 
tions, as the overall normalisations of F-^i and Flj>// are 
similar, meaning that the recoil spectrum of O 15 is ex¬ 
pected to differ strongly from the standard case in both 
energy and directional spectra. These two example distri¬ 
butions represent the most extreme departures from the 
standard SI/SD operator case and allow us to explore the 
full range of behaviours of the NREFT operators. 



^15,15 OC 


vlFs>+2- 


IV. STATISTICAL TESTS 


In this section, we explore how the two operators Oj 
and Oi 5 differ from the standard directional rate with 
regards to two statistical tests which have been proposed 
to confirm the DM nature of a directional signal. First, 
we explore how many events are required to reject the 
isotropy of the signal. Second, we determine how many 
events are required to confirm the median recoil direction 
of the signal. This allows us to quantify the particle 
physics uncertainties associated with the signals, arising 
from a lack of knowledge about which NREFT operators 
mediate the WIMP-nucleon interaction. 

We use the same experimental parameters as in 
Sec. |IIIB and assume perfect energy resolution and an¬ 
gular resolution. Of course, realistic experiments have 
finite energy resolution and are expected to have an an¬ 
gular resolution in the range 20 ° — 80°, depending on the 
recoil energy [60]. However, this idealised case allows us 
to place a lower limit on the number of events required to 
distinguish the DM signal from an isotropic background. 
Furthermore, our focus is on comparing different opera¬ 
tors, and assuming an idealised experiment allows us to 
disentangle experimental uncertainties from effects aris¬ 
ing from varying particle physics. 


A. Rejecting isotropy 

Backgrounds of terrestrial origin are expected to be 
isotropically distributed, so the first step in confirming 
the WIMP origin of a signal is to reject isotropy of the 
signal. We follow Morgan et al. [B2] and use the modi¬ 
fied Rayleigh-Watson statistic W*, defined in Ref. [HH] . 
as a measure of isotropy. Large values of W* indicate a 
larger degree of anisotropy. We calculate W* from the di¬ 
rections of mock events (distributed according to Eq. [^, 
discarding information about event energies. In order to 
determine the number of signal events iVwiMP required 
to reject isotropy, we use the following procedure for each 
value of VwiMp: 


FIG. 5. Number of WIMP signal events required to 
reject isotropy at 95% confidence. Results are shown 
assuming the signal is distribution according to each of three 
different NREFT operators: O4 (black •), O7 (blue T) and 
Ois (green ■). A Fluorine-based detector with an energy 
threshold of 20 keV is assumed. 


1. Generate 10000 mock data sets, each consisting 
of A^wimp recoil directions, distributed assuming 
a particular NREFT operator; 

2. Calculate W^jjyjp for each mock data set; 

3. Calculate the 5% percentile of 

4. Generate a further 10000 mock data sets, each 
consisting of Awimp recoil directions, distributed 
isotropically; 

5. Calculate for each mock data set; 

6 . Calculate the 95% percentile of 

The 95% percentile of >Vi*o is the value of W* above 
which we would reject isotropy at the 95% confidence 
level. We then find the value of Awimp for which this 
value is equal to the 5% percentile of W^jjyjp. For this 
value of Awimp, we can expect to reject isotropy at the 
95% level in 95% of experiments in which the signal 
events are distributed according to the NREFT opera¬ 
tor of interest. 

Figure]^ shows the results for the three operators O 4 
(the standard SD operator), O 7 and O 15 as a function 
of the WIMP mass. The number of events required to 
reject isotropy, assuming the standard SD operator, is 
of order 10 , in good agreement with the previous results 
of Morgan et al. [62]. As expected, O 15 requires fewer 
events to reject isotropy than the standard operators, 
while Or requires more. For example, at my, = 100 GeV, 
the standard operator O 4 requires 12 events to reject 
isotropy, compared to 16 and 11 events for Or and O 15 
respectively. Though this difference is relatively small 
in absolute terms, it represents an uncertainty of around 
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25% in the number of events required, arising entirely 
from particle physics uncertainties. 

At low masses, the number of events required for the 
three operators converges, with all three requiring only 
^ 6 events for a WIMP mass of 10 GeV. This is because 
at low WIMP masses, only large values of Wniin contribute 
to the event rate. This means that the exponential terms 
in Eqs.[^and[^will decay rapidly away from the forward 
direction, leading to highly directional rates for all three 
operators. However, as we increase the WIMP mass, the 
number of events required for each operator begins to 
diverge. In particular, A^wimp rises more rapidly for the 
operator O 7 . As demonstrated in Fig.[^ the maximum in 
the directional spectrum moves further from 0 = 0 as we 
increase leading to a more isotropic distribution of 
events and therefore more signal events required to reject 
isotropy. 


B. Confirming median recoil direction 

Once the isotropy of the signal has been confirmed, 
it will then be necessary to determine whether the me¬ 
dian recoil direction matches that expected from a WIMP 
signal. For all three operators we consider, the expected 
median recoil direction is the same and is in the direction 
of Flag. However, the distribution of observed median re¬ 
coil directions over an ensemble of experiments will be 
different. In order to quantify this, we follow Ref. [M] 
and examine the distribution of A, defined as the an¬ 
gle between the observed median recoil direction and the 
direction of Solar motion. 

As in Sec. |IV A[ we generate 10000 mock data sets for 
each hypothesised signal, as well as for the null hypoth¬ 
esis of isotropic recoils, and calculate the distribution of 
A. We then calculate the value of Awimp for which the 
5% percentile of A under the null hypothesis matches the 
95% percentile of A under the signal hypothesis. The re¬ 
sults are shown in Fig. 

As in the case of rejecting isotropy, the results agree 
with those presented previously |64j . with around 30 
events required to confirm the median recoil direction. 
Once again, O 4 and O 15 give almost identical results at 
low masses, as both have a highly directional recoil spec¬ 
trum. In this case, however, O 7 requires a larger number 
of events, approximately 50% more at low masses, in¬ 
creasing substantially as the WIMP mass is increased. 
Even at low mass, where the peak recoil direction co¬ 
incides with 0 = 0 , the directional spectrum is broad¬ 
ened by the structure of the transverse Radon Transform, 
arising from the coupling to This increases the size 
of fluctuations in the median recoil direction away from 
0 = 0. As the WIMP mass is increased, the ring-like 
structure described in Sec. |HI B] becomes significant and 
more events are required to confirm the median recoil di¬ 
rection. Above around 1000 GeV, the number of events 
required is around a factor of 2 higher than for the other 
two operators. 



FIG. 6. Number of WIMP signal events required to 
confirm median recoil direction at 95% confidence. 

Results are shown assuming the signal is distribution accord¬ 
ing to each of three different NREFT operators: O4 (black •), 
Or (blue T) and Oir, (green ■). A Fluorine-based detector 
with an energy threshold of 20 keV is assumed. 

While our results for standard interactions are in agree¬ 
ment with previous results, we have demonstrated that 
for NREFT operators, differences in directional spectra 
lead to different numbers of signal events required to con¬ 
firm the median recoil direction. For a 100 GeV WIMP, 
this required number of events ranges between 25 and 50, 
effectively introducing a factor of 2 particle physics un¬ 
certainty into otherwise model independent statements 
about the WIMP origin of a signal. 


V. COMPARISON WITH STANDARD 
INTERACTIONS 

Finally, we consider the possibility of discriminat¬ 
ing between NREFT operators and standard SI/SD op¬ 
erators. In Sec. EYi we considered relatively model- 
independent statistical tests which allowed us to distin¬ 
guish between signal and background. This was possible 
because the background hypothesis had a fixed form, that 
of isotropically distributed recoils. When comparing two 
signal hypotheses, however, this form is no longer fixed. 
This is because the signal rate for a given operator de¬ 
pends on the WIMP mass, which we assume is a 'priori 
unknown. The signal rate may also depend on other un¬ 
certainties, such as in the astrophysical distribution of 
WIMPs, or in detector performance. However, we ne¬ 
glect these uncertainties in the present study, focusing 
on the idealised, fixed-astrophysics case. 

In order to make statistically robust statements then, 
we must compare the observed distribution of events with 
that expected from each operator for all possible values 
of the WIMP mass. In order to account for this, and to 
make use of as much information as possible, we perform 
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a full likelihood-ratio analysis The null hypothesis 
Hq asserts that the signal is due entirely to standard SD 
interactions. The alternative hypothesis is that there 
is some contribution from another one of the NREFT 
operators (either Oj or O 15 , which we consider one at a 
time). The fraction of signal events due to one of these 
NREFT operators is denoted A, while the remaining frac¬ 
tion (1 — A) of events arises from the standard SD oper¬ 
ator. The likelihood L{mx,A) is then the probability of 
obtaining the observed event energies and directions for 
a given value of and A. 

For each of the 10000 mock data sets, we generate 
.^wiMP events, distributed in energy and direction ac¬ 
cording to either Oj or O 15 . We then calculate the fol¬ 
lowing test statistic: 



go = -2 In 


L(m^,0) 

L(m^,A) 


(29) 


Here, L{m^, 0) is the likelihood under the null hypothesis 
(i.e. SD-only events), maximised over all values of m^. 
The unconditional maximum likelihood is then denoted 
L{rh^, A), where we maximise over both the WIMP mass 
and the non-standard operator fraction A. Accord¬ 
ing to Wilks’ theorem [ 66 ], for data distributed under 
the null hypothesis (SD interactions-only), go is asymp¬ 
totically x^-distributed with one degree of freedom (the 
difference in dimensionality between the null and alter¬ 
native hypotheses)]^ 

For each value of IVwimp, we calculate g®®^, defined 
such that 95% of experiments observe a value of go 
greater than or equal to go^^°. We can then calculate 
the p-value (and corresponding confidence level) for the 
SD-only hypothesis, based on this value, as 

noo 

P= -P(x?)dx?, (30) 

where P{xi) is the probability density function for the 
Xi distribution. That is, we calculate the probability of 
observing a value of go as large or larger than g®^^ if all 
signal events are due to standard SD interactions. A p- 
value this small or smaller will then be obtained in 95% 
of experiments in which the signal is due to non-standard 
interactions. The smaller the value of p, the greater the 
confidence level with which we can reject the SD-only 
hypothesis and infer the presence of other NREFT oper¬ 
ators. 

The resulting p-values and confidence levels (in units of 
a) are shown in Fig.j^as a function of A^wimpj assuming 
signals arising from O 7 (blue) and O 15 (green), with a 
WIMP mass of 100 GeV. We show the results obtained 
when the full energy and directional information is used 


® We have verified this numerically for several values of A'^vvimp 
above 30. 


FIG. 7. Confidence levels for rejecting standard SD- 
only interactions. We show the p-value obtained in 95% 
of experiments as a function of the number of signal events 
Awimp for the null hypothesis of standard SD-only interac¬ 
tions. We also show the corresponding confidence level with 
which standard interactions can be rejected (in units of a). 
The signal is distributed according to the operators O 7 (blue 
T) or £>15 (green ■), with a WIMP mass of 100 GeV. A 
Fluorine-based detector with an energy threshold of 20 keV 
is assumed. Experiments with and without directional infor¬ 
mation are shown as solid and dashed lines respectively. Note 
that the curve for O 7 using only non-directional information 
lies above p > 0.9 for all values of Awimp considered. 


to calculate the likelihood (solid lines) as well as those 
obtained using energy information only (dashed lines). 

In the case of non-directional detection, the p-value ob¬ 
tained for O 7 remains large (p > 0.9), even with up to 
1000 signal events. Figurej^shows the energy spectra for 
the operators considered in this section and illustrates 
that the spectra for the standard SD (O4) and O7 oper¬ 
ators are almost identical. As described in Sec. [Ill B[ for 
light nuclei such as Fluorine, we expect differences be¬ 
tween the different form factors to be negligible. The 
most significant difference between the two spectra is 
therefore that the velocity integral for O7 is p^(vniin), 
defined in Eq. 22 which is weighted by v]_. 

The behaviour of r]"^{vynin) can be understood by ex¬ 
amining Eq. 19 If cr„ is large compared to uiag, then the 
exponential term in the TRT will vary slowly as a func¬ 
tion of uiag • g = uiag cos 0. When we integrate over all an¬ 
gles, the term proportional to uiag -g is then sub-dominant 
compared to the remaining terms, meaning that the ve¬ 
locity integral will be roughly proportional to 

the standard integral r]{vjnin)- If instead is small, the 
exponential will be sharply peaked and the angular inte¬ 
gral will be dominated by iTiag • g = Vmin- In this case, 
the term proportional to Flag • 9 can be comparable to 
the remaining terms and will give an extra contribution 
to the energy spectrum proportional to —u^i„. 

Physically, a small value of ay leads to a narrow veloc¬ 
ity distribution. This means that almost all WIMPs are 
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FIG. 8. Energy spectra, normalised to unity, for sev¬ 
eral NREFT operators considered in Sec.[^ On (black 
•), O 7 (blue T) and O 15 (green ■). The energy spectra are 
normalised to a single event in the range Er £ [20, 50] keV. 
We assume = 100 GeV and a Fluorine detector. 


travelling with velocity close to wiag. Due to the coupling 
to scattering in the forward direction is suppressed, 
meaning that scattering through larger angles (and there¬ 
fore with lower recoil energies) is enhanced. This is seen 
in the right panel of Fig. and was previously discussed 
in Sec. |III A[ For larger values of cr„, there is a significant 
population of WIMPs travelling with large velocities at 
an angle to Flag. These can induce high energy recoils in 
the direction of Flag while still satisfying the kinematic 
constraint v ■ q = Wmin, meaning that the energy spec¬ 
trum at high Ei 7 is not depleted. This can be seen in 
particular in the forward direction in the right panel of 
Fig. 0 In the case of the SHM, with uiag = y/2o-v, this 
effect means that the energy spectrum obtained from the 
TRT matches the standard case within a few percent. 
Such small differences can easily be compensated for by 
varying the WIMP mass, making discrimination difficult 
with energy-only detectors. 

When directional information is included, however, the 
distribution of events can be distinguished. As is clear 
from Fig. the transverse Radon Transform leads to 
a different angular distribution of events for Oy when 
compared with the standard RT for the SD interaction. 
These distributions are sufficiently different that with 
around 500 events, the standard interactions can be re¬ 
jected at the 2(7 level in 95% of experiments which are 
directionally sensitive, with 3a discrimination possible 
with around 700 events. We emphasise that in this case, 
energy-only information does not allow us to significantly 
distinguish between the two operators. Thus, it is only 
the directionality of the signal which allows us to dis¬ 
criminate and reject the SD-only hypothesis, in favour of 
the NREFT operator O 7 . 

By contrast, energy-only experiments can very quickly 
distinguish a signal dominated by O 15 from a standard 


SD-only signal. The 2a level is reached in 95% of exper¬ 
iments with around 150 signal events, while the 5a dis¬ 
covery level could be achieved with as few as 550 events. 
This is because of the characteristic energy spectrum pro¬ 
duced by Oi 5 , which rises as g® for small g, before form 
factor suppression becomes important and the spectrum 
flattens at high energy. This spectrum is also illustrated 
in Fig. and cannot be easily mimicked by the standard 
SD signal, even if the WIMP mass is varied. 

When directional information is also included, again 
we see a significant improvement in the confidence level 
with which standard interactions can be rejected. For 
a given number of signal events, the rejection of stan¬ 
dard interactions is approximately Icr more significant 
when directional information is included. As a result, a 
5(7 rejection of SD-only scattering can be achieved for 
around 300 signal events. This arises because the nu¬ 
clear response function for O 15 , shown in Eq. 27 con¬ 


tains a contribution proportional to q'*v]_ as well as a 
contribution which goes as g®. This means that as well 
as producing a different energy spectrum, O 15 also pro¬ 
duces a different directional distribution. Because the 
two contributions to the nuclear response function have 
a similar normalisation, this difference in directionality 
can be easily observed and can be used to distinguish 
Oi 5 from the standard SD-only case. 


VI. DISCUSSION 


In this work, we have considered the directional recoil 
spectra produced by the set of NREFT operators. We 
have focused on a single, idealised Fluorine-target detec¬ 
tor for concreteness. Though we have assumed a reason¬ 
able energy threshold of 20 keV [59], we have assumed 
perfect angular resolution (compared to the typical reso¬ 
lution of 20°-80° jSD]). We have also assumed that there 
is no background contamination. Introducing a finite an¬ 
gular resolution and background would increase the num¬ 
ber of events required to reject isotropy and confirm the 
median recoil direction. The numbers reported above in 
Sec. IV therefore represent a lower limit and illustrate 
that even with an idealised detector, uncertainties com¬ 
ing from particle physics can be as much as a factor of 
2 . 


As well as experimental uncertainties, these results are 
subject to astrophysical uncertainties. Though we have 
presented the TRT for a stream distribution for illustra¬ 
tion purposes in Fig. we have restricted calculations 
to the SHM with fixed parameters. In Ref. [H2], the im¬ 
pact of different halo models on the number of events 
required to reject isotropy was studied. The values ob¬ 
tained therein vary by around 20 %, meaning that the 
particle physics uncertainties presented here are expected 
to be comparable to astrophysical uncertainties. 

We note in particular that one of the conclusions of 
SecjV]- that O 7 and the standard SD operator O 4 are in¬ 
distinguishable without directional information - is only 
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true for SHM-like velocity distributions. In cases where 
Gy is significantly smaller than uiag (such as in the pres¬ 
ence of a stream), the recoil spectra for the two operators 
will diverge. However, this is also likely to accentuate the 
differences between the angular event distributions for 
the two operators, meaning that we expect that direc¬ 
tional sensitivity will still provide a significant improve¬ 
ment in discrimination. The impact of astrophysical un¬ 
certainties in the NREFT framework has been briefly dis¬ 
cussed in the past [3T], though clearly a detailed study 
of such uncertainties in directional and non-directional 
experiments will be necessary in future. 

Even if the velocity integrals for the operators O 4 and 
Oj are indistinguishable, discrimination between these 
different operators may also be possible through other 
methods. For heavier nuclei, such as Xenon, there may be 
more significant differences in the form factors associated 
with each operator. By comparing the energy spectra 
and number of events in several experiments using differ¬ 
ent target nuclei, it may be possible to determine which 
form factor (and therefore which operator) is mediating 
the interaction (see e.g. [2S])- However, many NREFT 
operators lead to interactions with the same form factor. 
In addition, uncertainties in calculations of the form fac¬ 
tors and in the value of the WIMP mass may make such 
an approach more difficult EIIIT]. A more promising ap¬ 
proach is to measure the annual modulation of the dark 
matter signal, which has a different time-dependence for 
operators coupling to w_l. However, the annual modula¬ 
tion would have to measured in several experiments and 
compared before different operators could be discrimi¬ 
nated [ 68 j . The use of directional information instead 
allows operators to be distinguished with just a single 
experiment. 

Finally, we caution that the likelihood-based approach 
of Sec. [V| would not be appropriate for all NREFT oper¬ 
ators. As previously discussed in Sec. |HIB| those opera¬ 
tors which differ from the standard SI/SD interactions 
only through suppression do not change the direc¬ 
tional dependence of the full double-differential recoil dis¬ 
tribution (Eq. 151, but only affect the energy spectrum 
of events. Thus, the addition of directional information 
does not improve prospects for discrimination compared 
to the energy-only case, as the directional dependence of 
both operators is the same. In spite of this, the statisti¬ 
cal tests of Sec. m are still useful for such operators, as 
they allow robust and relatively model-independent com¬ 
parisons to be made with isotropic backgrounds, without 
any reference to the energy spectrum of the operator. 


VII. CONCLUSIONS 

In the current work, we have explored the directional 
signatures of DM-nucleon interactions within the frame¬ 
work of non-relativistic effective field theory (NREFT). 
Some of the operators arising in NREFT lead to a sup¬ 
pression by the recoil momentum and therefore can 


affect the angular distribution of recoils when the di¬ 
rections (but not energies) of the events are considered. 
Other operators lead to a coupling to the WIMP velocity 
perpendicular to the recoil direction v\. In these cases, 
the directionality of recoils can be fundamentally differ¬ 
ent compared to the standard SI and SD interactions. A 
number of the NREFT operators give similar directional 
recoil spectra and can be classified, as in Eq. accord¬ 
ing to the scaling of their cross sections with and 
We have focused on two operators in particular, 
namely 


Or = Sn ■ V , and 


O 


15 


= -(4-:f)((*5nXU^).^). 


(31) 

(32) 


The operator O 15 produces recoils which are more 
strongly peaked in the forward direction, compared to 
the SI/SD case. The operator Or instead produces recoils 
which are the most isotropic of all the NREFT operators. 
In Eq. |17[ we have defined the transverse Radon Trans¬ 
form (TRT), which takes into account the u^^-weighting 
of the interaction cross section, and which appears in the 
recoil spectrum for Or as well as several other opera¬ 
tors. The TRT suppresses scattering in the forward di¬ 
rection and can therefore produce a pronounced ring-like 
feature in the directional recoil spectrum. Though ring¬ 
like features have previously been discussed in the con¬ 
text of standard operators m, the ring arising from the 
TRT should be observable down to lower WIMP masses 
(tojj. < 20 GeV) and for higher threshold energies (as 
high as 20-30 keV for high mass WIMPs). 

In Sec. IV we have explored the number of events re¬ 
quired for each operator to reject the isotropy of the re¬ 
coils and to confirm the median recoil direction. For the 
strongly directional operator O 15 , only a slightly smaller 
number of recoils is required compared with standard SD 
interactions. For Or, however, substantially more recoils 
are required, due to the broader recoil distribution in this 
case, with the difference increasing as we consider higher 
WIMP masses. For a 100 GeV WIMP, we conclude that 
10-15 events are required to reject isotropy, while 25-50 
events are required to confirm the median recoil direc¬ 
tion, depending on the operator in question. 

In Sec. |Vj we considered how well the NREFT op¬ 
erators could be distinguished from the standard case 
using both directional and non-directional information. 
For an underlying O 15 signal, SD-only interactions could 
be rejected at the ia level with around 300 events if 
only energy information is available. Including informa¬ 
tion about the recoil directions, this number is reduced 
to around 150 events. For Or, the energy spectrum of 
events cannot be distinguished from the standard case. 
However, when differences in the angular distribution of 
events are taken into account, it may be possible to ex¬ 
clude the SD-only scenario at the 3 ct level with around 
700 events. Though such large numbers of events would 
require large exposures (and are certainly beyond the 
scope of current experiments), future experiments with 
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larger target masses and lower energy thresholds may be 
able to distinguish the different operators. 

We have demonstrated that directional information 
may be the only means of discriminating those operators 
which couple to iTj_ from those which do not. Though 
this study is far from exhaustive (neglecting, for exam¬ 
ple, interference terms between different operators), we 
have highlighted the importance of directional detection 
for probing the particle physics nature of Dark Matter. 

During the preparation of this manuscript, a pre-print 
also discussing the directional rates in NREFT was made 
available online (Ref. [SS] ). In that paper, the author con¬ 
siders the directional spectra and relative contributions 
of different NREFT operators for several possible tar¬ 
get materials in directional detectors. Instead, we have 
considered a single target material (CF 4 ) and focused 
on comparing the directional spectra produced by each 
operator. We have also considered the possibility of dis¬ 
tinguishing between different operators using directional 
detection. However, the results of this paper and Ref. [51] 
are in broad agreement, including our expressions for the 
Radon Transforms and our predictions of a novel ring-like 
signature for certain NREFT operators. 
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Appendix A: Non-relativistic operators 

Here, we list the NREFT operators which are consid¬ 
ered in this work. At the nucleon level, they are con¬ 
structed from the following Hermitian operators: the mo¬ 
mentum transfer the transverse WIMP-nucleon 

velocity u-*-; the DM spin and the nucleon spin Sn- 
The transverse velocity is given by 


vj_ = V -\—— , (Al) 

P-X" 

where /i^„ = m^mn/(m^ -\- irin) is the WIMP-nucleon 
reduced mass and is the nucleon mass. The list of 


NREFT operators is then as follows [221HH [5^ : 

= 1 

O3 = iSn ■ (— X V^) 

rrin 

04 = 4 • 4 

O5 = iS^ ■ (-^ X ff-*-) 

Oe = (4 ■ 91(4 • q) 

O 7 = 4 • 

Os = Sx- 

Og = iS^ ■ (4 X q) (A 2 ) 

Oio =iSn- q 

Oil = iSx ■ q 

O 12 = 4 • (4 X v^) 

Oi3 = i(4-^^)(4-—) 

TTln 

Oi4 = i( 4 -—)(4-^7^) 

TRn 

Oi5 = -(4 -—)((4xl;^)-^). 

We neglect the operator O 2 = 4) it does not arise 
at leading order from a relativistic Lagrangian without 
significant cancellation. It is therefore typically sub¬ 
dominant to other operators in the list. We also omit 
the two operators recently reported in Ref. |2S|. Several 
dictionaries are available which allow one to translate 
from a relativistic interaction Lagrangian to the NREFT 
operators listed above |23| |25l |56| . 

In addition, we have considered an example of a long- 
range operator: 


= ^ . (A3) 


This operator behaves as Oi, with an additional q sup¬ 
pression of the nuclear response function. 

From the nucleon-level operators, it is necessary to cal¬ 
culate the matrix elements of these operators within the 
nucleus, summing over the contributions of all nucleons. 
Neglecting interference terms between different opera¬ 
tors, the resulting nuclear response functions are given 
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by: 

Fia = Fm ■ 


-^ 3.3 = ( v\F^' + 


Faa = 


-^5,5 = 


1 

8 m' 

g(jx) 

16 
CiJx) 


m- 

(_Fe/ + Fy,") , 


T 


v_\_Fm H-5"-Fa 


p C'Ox) p 

-F"?,? = g^^i-F"i;/ , 
C(jx) 


-Fg.s = 


v]_Fm + -^-Fa 




P ^ C(jx) g" p 

-^9,9 -| ^ 9 -^Yl' 

lo 

-f^lO.lO — 7—7-f^E" , 

4 


1 

F"ll,ll = 7— ^Fm : 

4 m„ 


F"i 2 a 2 = 

-^13.13 = 
Fi4,14 = 
-F’ls.is = 


Cjjx) 

16 

CjJx) 

16 

C(jx) 

32 m2 

g(jx) g" 

32 m^ 


F’e" + + —7 (F^$/ + Fqyfi) 


v\F^ii + ^F^, 
mi 


v^Fy.' 


v±Fy' + 2 — 7 -F$' 
mi 


(A4) 


Here, C(j;^) = 4j^(j^ + l)/3, where is the DM spin. 
The transverse velocity appearing here is the (complex- 
valued) WIMP-nucleus velocity: 

vj_ =v-\ - — , (A5) 

^J^xN 


with N denoting the nuclear^ rather than nucleon^ mass. 
We have suppressed the isospin indices (r, r') which ap¬ 
pear in Eq. [T^ 


The functions Fm, Fy', Fy", Fa, F^, and F$" are the 
standard nuclear form factors appearing in the study of 
semi-leptonic electroweak interactions [70] and are func¬ 
tions only of q^ for a given nucleus. For the ^^F nucleus, 
we use the form factors given in Ref. [23] obtained un¬ 
der the one-body interaction approximation. Under this 
assumption, the form factors decay approximately expo¬ 
nentially with ^2, leading to suppression at large recoil 
energies. 

We assume isospin-zero interactions in this work, so 
only those form factors with r = r' = 0 will be relevant. 
In order to compare the relative strengths of the different 
form factors, we report below their values for at g = 
0 : 


Fm(0) = 90.25 
Fy>{0) = 0.435 
Fs-(O) = 0.218 
Fi,i0) = 0 
F$-(0) = 0.123 
Fa(0) = 0.0015 


(A6) 


This means that where a response function contains con¬ 
tributions from two or more form factors, some may be 
subdominant. In the case of Og, for example, there are 
two terms in the recoil spectrum - one coupling to v'^ 
and the other to q'^. However, the form factors associ¬ 
ated with each term {Fm and Fa) differ in normalisation 
by roughly 5 orders of magnitude. This means that the 
q2-term in Fg^g can effectively be neglected. 
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